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An approximate method for determining the effect of the levelj 
of turbulence on the aerodynamic characteristics of con- 
vergent and diffuser channels is examined. Utilized for thi^ 
purpose is a pulse equation for the boundary layer, in which 
the turbulence of the external flow is introduced on the 
basis of experimental data, which represent the coefficient 
of friction c f and the form-parameter Ho on the degree of 
turbulence E. It is found that, at significant levels of 
external turbulence, losses must be considered not only in 
the boundary layer, but also in the central region of the 
channel . 
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LOSSES IN CHANNELS WITH INCREASED EXTERNAL TURBULENCE 


A. Ye. Zaryankin and G. S. Solov'eva 


Examined herein is an approximate method for calcu- /164 * 

lating the effect of the level of turbulence on the aero- 
dynamic characteristics of convergent and diffuser chan- 
nels. Utilized for this purpose is a pulse equation for 
the boundary layer, into which the turbulence of the ex- 
ternal flow is introduced, based on experimental data 
which represent the coefficient of friction c^. and the 
form-parameter Ho on the degree of turbulence r E. 

The effect of external turbulence on the aerodynamic charac- 
teristics of channels has been studied little, although, in a 
number of cases, the effect of the turbulence of an approach stream 
has proved decisive. 

In particular, by examining the determination of the losses 
in turbine cascades, pipes and diffusers, it is entirely necessary 
to introduce the degree of turbulence of the external flow into 
the calculation, and to consider this magnitude one of the de- 
cisive parametric criteria of similitude in the experiment. 

The following considerations graphically emphasize what has 
been said. Let a flow with an increased degree of turbulence E 
and initial parameters p 0 and to (Fig. l,a), with a pressure be- 
hind the cascade p 2 , flow into a turbine cascade. We will select 
a current line Oa in the channel between the vanes, for which the 

local coefficient of losses t. will have a minimal value. For the 

1 

other current lines Ob, Oc, and so on, ? > ? i min . The process of 
expansion in the thermal diagram, which corresponds to the selected 
lines Oa, Ob, Oc, is presented in Figure l,b. By expressing the 
local loss coefficient in standard form, in parts of the avail- 
able energy in the outlet section, we obtain 

£i = l-ff,/ffo = l-C 2 i 2 /C2! 2 , 


*Numbers in the margin indicate pagination in the foreign text. 
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Fig. 1. Flow process in the cascade with increased 
turbulence of the external flow (a) and its de- 
piction in the is-diagram (b) . 


here, c 2i is the outlet velocity of the flow at the examined point, 
and c 2t is the theoretical velocity, determined by the isoentropic 
differential H 0 . 


We will average the magnitudes of according to mass 


PltCltFl f P2i<?2< / t - 2( 2 \ F 

fcs I (1 | d . 

m J P21C21 \ c u 2 } Fz 
0 


( 1 ) 


The product P2t C 2t F 2 =m t corres P oni ^ s to the theoretical flow rate /165 
with the absence of losses. We will use Am to designate the diff- 
erence in the flow rates 


A m = mt — m = 


1 

_ C I . p 2 i c 2 t \ _ 

p2(C2<^2 I I 1 | dF. 

“ \ P2(C2i / 


( 2 ) 


Then 

m=m(-Am=p2(C2(/ ? 2 [ 1 - f / 1 - — \ 1 dF. 

L J \ pztCit / J 

The integrals in relationships (1) and (2) coincide in structure with 
those expressed for a relative area of energy losses <$ 2 ***, and a 
relative area of displacement <$ 2 *' hut the dimensionless velocities 
c 2 - are expressed in parts of the theoretical velocity c 2t , rather 
than parts of the velocity c 2i in the examined section. 


By analogy with the magnitudes 6* and 6***, we will consider 
the integrals 
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A^r/i w (3a) : 

J \ P2(C2( / ( 

'° ! 
... f P atC2< / 1 \ 

J p2(C2( \ C2( 2 / 


dF ( 36 ) j 


as the generalized areas of displace- 
ment and energy loss. 

Actually, with the absence of 

losses on the zero current line Oa, 

c. =c„ and A - *=6, while ~K * * *= 

_2t 2max 2 

5 * * * 


If one introduces the co- 
efficient of velocity for the zero 

current line <P<>, equal to <po=c„ /c 

- 2max 



8 to 


Fig. 2. Dependence of co- 
efficient of losses on the 
coefficient of velocity 
at the zero current line. 
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into the examination, then one may 

establish the following simple relationships between the generalized 
and standard areas 62 * and 

A2*— 1— <po (1- 62*) ; A2 “(po(l — 62*) — 

-<Po , (l-6l'-S>')- (4) 


By substituting expression (2) into formula (1), with regard for 

the introduced designations (3) and (4), we obtain 

•• 

r 62 i 

?e — 1 — <Po I 1 I • (5) 

L i-6 2 ' J 

With E=0 , q?o =1 and 

Thus , 

U = l-tpo(l-U-o). (6) 

Relationships (5) and ( 6 ) are common for evaluation of the losses 
in random channels. We would note that, during their derivation, 
we put nature aside in decreasing the coefficient of velocity 
on the line of minimal losses. This may be increased turbulence 
of the approach stream, non-uniformity, periodic nonstationariness , 
and so on. 
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How strongly the coefficient of losses varies, as a function 
of the causes which evoke a reduction in q>o. , is evident from Figure 
2, where the coefficients of losses £ =0, with the presence of a 

C 

potential focus, are plotted along the x-axis and the value of the 
losses which corresponds to the coefficient <P<> . if one assumes 
that the integral areas of the boundary layer do not depend on E, 
then, even with <p 0 (=0.98, the losses differ substantially from the 
losses found with a zero or small degree of turbulence of the 
approach stream. 

In actuality, A 2 * and & 2 *** vary appreciably with a change in 
turbulence, or its increase leads to an increase in the completeness 
of the profile, the Reynolds stresses, and an increase in the physi- 
cal thickness of the boundary layer. 

Thus, the problem of evaluating the losses in the channels 
amounts to the determination of the coefficient <Po,, which charac- 
terizes the dissociation of energy in the nucleus, and the calcu- 
lation of the minimum areas of the boundary layer, which develop /166 
under conditions of increased turbulence in this nucleus. 

The question of evaluation of the coefficient cp 0 presents 
an independent problem, the solution of which is the subject of 
studies [1,2] . 

Specifically, based on the statistical theory of turbulence 
in [1] , it is shown that, for a pipe, the magnitude of <Po may be 
found according to the relationship 

<p 0 =yi— 9,55E 3 T. 

The experimental determination of this coefficient in turbine 
cascades is given in [2]. 

The second part of the problem, which consists of evaluating 
the magnitudes of and ~§*** , is more complex and, in this stage, 
may be solved with a great deal of approximation. One of the po- 
tential means of this evaluation comes from the properties of the 
solution of Carman's integral relationship. 
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For simplicity, we will examine the case of a flat boundary 
layer, and we will limit ourselves to a comparatively low level of 
turbulence (E<7-8%), when one may still disregard the additional 
terms in the integral equation, and utilize its common form for 
analysis 


d6“ u t '6" t „ c,- 

+ (H+ 2 )= , 

dx «i pUi 2 2 



where Ui is the velocity of the external flow and ui’=dui/dx. 

We will introduce the dimensionless magnitudes c^ and H 

(c =c~ c,., H=H 0 H) into the examination, taking the magnitudes c- 
± t o i ±o 

and H 0 as the norming factors, which correspond to a gradientless 
flow, for which Ho=1.4 and 

c/ 0 /2=£o/Re*’ m . 


With regard for what has been said, equation (7) takes on the form 

<Z6” ui'6** 

Re" m + Re ,,m (2+// 0 5') = c,l o. (8> 

dx Ui 

Simple transformations lead to the relationship 
d 

— (Re-^'^tTO+lKoc/ +[m-( m+ i) (2+HqH) ]T, 
dx 

where 

ui'6** 

r Re" m . 

Ui 

In examining the question of the change in the normed magni- 
tudes of c f and H, one may note that c f hardly changes at first, 
according to the available test data in the diffuser region, i . e . , 
it remains close to one, and then drops sharply to zero at the 
point of break-away of the boundary layer. In turn, the parameter 
H increases continuously right up to break-away of the flow. 

By utilizing the nature of the change in the indicated magni- 
tudes from the complex r, we will approximate and H with the 
following relationships: 

c/=l-fc,r 2 ; H=>l-k 2 T+k 3 T 2 . 
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We will show that the accuracy of the proposed approximation 
(numerical values of the magnitudes ki , k 2/ k 3 ) does not affect the 
quadrature of Carman's equation. Actually, by substituting (10) 
into (9) , we obtain 
d 

(6** Re” m ) — a + bT + eT 2 + dT 3 . (11> 

dx 

Here, a=(m+l)c 0 ; b=- [ (m+1 ) (1+H' 0 ) (1+H )+l]; e= (m+1 ) (k 2 H 0 -C ok ); d= 

(m+1 ) H 0 k 3 . 

The polynomial in the righthand part of equation (10) differs 
very little from the linear function, since the coefficients e and d, 
which are small per se , stay factors at high levels of the parameter 
r, the magnitude of which does not exceed 0.04-0.06 at the break- 
away point. In other words, without committing a large error, one 
may linearize relationship (11) and, as a result, obtain the well- 
known quadrature 

(12) 

are equal to /167 

a 2 = m </(m + 1); 

a t =l/(m + l). 

The obtained conclusion makes it possible to reveal an impor- 
tant property of the solution of (12), the essence of which consists 
of the fact that all of the constant coefficients in (12) are com- 
pletely determined by the magnitudes which characterize the flow on 
a plate. 

From here, it follows that the accuracy of calculation of the 

boundary layer in a gradient flow, based on the equation of T. Carman, 

depends on how accurately the dependences c_ and Ho are determined 

f o 

for a plate 1 . 


1 This method is utilized, in essence, in the known layer-calculating 

method of K. K. Fedyaevskiy, A. V. Kolesnikov and A. N. Smolyaninova , 

where the increased accuracy is achieved because of the introduction 
of more specific data on a gradientless flow [7] . 


&** * 

0 a x Pf ^ % 1 a A 

L = R e?'u2«+ a .«< |J ■ 

x 

The coefficients of equation (12) 

*J3=1 + (to + 1) (/To+1) » 
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By utilizing this very method and the examined property of 
equation (1) , it is comparatively easy to introduce the level of 
turbulence of the external flow into the calculation. For this 
purpose, it is sufficient, as has already been noted, to carry out 
the detailed study of the boundary layer on a plate with an in- 
creased value of the magnitude E, and to utilize the obtained data 

on the coefficients c~ and Ho for calculation of the values of 

t o 

in equation (12) . 



E.% 

Fig. 3. Dependence of the magnitude of 6 E **/ 

6** 0 on the level of turbulence with convergent 

(u=0.3 + 0.7x) (1) and diffuser (u=l-0.7x) (2) flows. 


The still evidently insufficient available experimental data 

on the dependences of c £ and Ho on the level of turbulence E [3-6] 

t o 

may be roughly approximated in the following manner: 


c,*=c /o ( 1+0, 01£ 2 ), 

/ 0,25 \ 

Ho B = Ho ( 0,75 + — ) 

V 1+0,1 E) 


(13) 


where E is the level of turbulence in percent. 


With regard for the listed expressions for the coefficients 
ai and a 3 , we obtain the following dependences: 


a, = 1 + (m 


ai — [(to + 1)£o(1+0,01£ 2 )], 

0,25 


(14) 


[ «. ( 0.7 5 + ) + '] 
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Their substitution into formula (12) gives 


6 =• 


("» + i) ; 0 (i +o,oi£ 2 ) 


1 1 , 2+™ . „ f 0,25 \ 

l m+1 (u!!M Tn + 1 ^" m +* ' ' ’ + l+o.i E> 


(Re^) m+ 1 (u !?) 1 


fi 


* 2+m+(m+1) H . ( 0 , 75 + 

\ui 


dx 


1 

m 4-i 


(15) 


j 


! 


The results of the calculation, according to relationship (15), 
of the integral depth of the loss of the pulse <$£* in the outlet 
section of the convergent and diffuser channels are represented as 
a function of the level of turbulence in Figure 3. 


We would note that the possibilities of the examined procedure 
are not limited to the case of increased turbulence of the external 
flow, but are considerably broader. Specifically, the very same 
method may be utilized to solve the problem of calculating the 
turbulent boundary layer in a flow of compressible gas. What is 
more, after the accumulation of data on the resistance of a plate 
and the integral characteristics of the layer by analogous means, 
one may correct the solution for calculation of moisture during the 
movement of a moist vapor flow in the cascades of turbines. 


Conclusions 

1. With considerable external turbulence of the flow, it 

is not necessary to take into account the losses, not only within 
the limits of the boundary layer, but also in the central part of 
the channel. 

2. Based on the integral relationship of Carman, we obtained 
an analytical expression for calculating the depth of displacement 
with an increased external turbulence of the flow. 


/168 
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